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The rich phase diagrams of magnetically frustrated pyrochlores have maintained a high level of 
interest over the past 20 years. To experimentally explore these phase diagrams requires a means of 
tuning the relevant interactions. One approach to achieve this is chemical pressure, that is, varying 
the size of the non-magnetic cation. Here, we report on a new family of lead-based pyrochlores 
A 2 Pb 2 0 7 (A = Pr, Nd, Gd), which we have characterized with magnetic susceptibility and specific 
heat. Lead is the largest known possible B-site cation for the pyrochlore lattice. Thus, these mate¬ 
rials significantly expand the phase space of the frustrated pyrochlores. Pr 2 Pb 2 0 7 has an absence 
of long-range magnetic order down to 400 mK and a spin ice-like heat capacity anomaly at 1.2 K. 
Thus, Pr 2 Pb 2 0 7 is a candidate for a quantum spin ice state, despite weaker exchange. Nd 2 Pb 2 0 7 
transitions to a magnetically ordered state at 0.41 K. The Weiss temperature for Nd 2 Pb 2 0 7 is 
dew = —0.06 K, indicating close competition between ferromagnetic and antiferromagnetic interac¬ 
tions. Gd 2 Pb 2 0 7 is a Heisenberg antiferromagnet that transitions to long-range magnetic order at 
0.81 K, in spite of significant site mixing. Below its ordering transition, we find a T 3 / 2 heat capac¬ 
ity dependence in Gd 2 Pb 2 0 7 , indication of a ground state that is distinct from other gadolinium 
pyrochlores. These lead-based pyrochlores provide insight into the effects of weakened exchange on 
highly frustrated lattices and represent further realizations of several exotic magnetic ground states 
which can test theoretical models. 


INTRODUCTION 


Magnetic frustration occurs when it is impossible to 
satisfy all pairwise magnetic interactions simultaneously 
and is typically a consequence of the underlying lattice 
geometry. The canonical example of magnetic frustration 
is antiferromagnetic Ising spins placed on the vertices of 
an equilateral triangle; in this scenario, only two of the 
three spin pairs can be anti-parallel, resulting in mag¬ 
netic frustration. Most structures studied in the context 
of magnetic frustration are composed of triangular mo¬ 
tifs [lj. Some examples of frustrated lattice geometries 
include the kagome, honeycomb, and perovskite lattices. 
However, there is perhaps no lattice better suited to the 
study of magnetic frustration than the pyrochlore lattice 
[2]. Pyrochlores, with formula A 2 B 2 O 7 , are composed of 
two interpenetrating networks of corner-sharing tetrahe- 
dra. The vertices of one tetrahedral network are occupied 
by the A 3+ cations while the second network is occupied 
by the B 4+ cations. This corner-sharing tetrahedral ar¬ 
rangement is subject to extreme frustration when occu¬ 
pied by a magnetic ion. It is fortunate for the study of 
magnetic frustration that a large fraction of the cations 
which can occupy the pyrochlore lattice are, in fact, mag¬ 
netic. 

The most well-characterized family of frustrated 
pyrochlores are the rare earth titanate pyrochlores, 


R 2 Ti 2 0 7 [U- The B-site cation, Ti 4 +, has a closed va¬ 
lence shell and is thus non-magnetic. By varying the 
A-site rare earth cation, many unique and exotic mag¬ 
netic ground states have been uncovered in this family. 
The spin ice state, a magnetic analog for the proton ar¬ 
rangement in water ice, is observed for R = Ho 0, Dy 
0. WWhen the rare earth site is occupied by R = Tb, 
there is an absence of conventional long-range magnetic 
order down to 70 mK, leading Tb 2 Ti 2 0 7 to be considered 
a spin liquid candidate 0, Q ■ Long range order may also 
be absent in the case of the apparent quantum spin liquid 
R = Yb 0 , 0 . A detailed investigation of the Hamilto¬ 
nian of Yb 2 Ti 2 0 7 has justified the label “quantum spin 
ice” [§]. However, the low temperature magnetism of 
Yb 2 Tj 20 7 remains controversial, due to sample variation 


10, 11] and conflicting reports [12|, 13]. Even the ti¬ 


tanate pyrochlores that transition to long-range magnetic 
order, R = Gd, Er, have exotic properties. In the case of 
Er 2 Ti 2 0 7 , one proposed scenario is that the transition to 
long-range magnetic order is selected by thermal fluctu¬ 
ations (i.e., order-by-disorder) GHH. I n Gd 2 Ti 2 0 7 , an 
uncommon example of Heisenberg spins on a pyrochlore 
lattice B , the antiferromagnetic ground state is a com¬ 
plex multi-fc structure with a two-stage ordering BO. 


The plethora of exotic magnetism discovered in the 
titanate pyrochlores was further expanded through sub¬ 
stitution of titanium for other non-magnetic cations of 
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varying size, i.e., chemical pressure [20j, 21]. Studies 


with chemical pressure have allowed investigation into 
the phase space for these frustrated materials. In the py- 
rochlores, variation of the non-magnetic B-site cation has 
two main consequences. Firstly, the change in ionic radii 
introduces a positive or negative chemical pressure. This 
results in differences in the distances between magnetic A- 
site cations, allowing an effective tuning of the exchange 
interaction [H]. The second consequence of substituting 
the B-site cation involves subtle changes to the local oxy¬ 
gen environment surrounding each magnetic cation [231 ]. 
These differences, although structurally small, can result 
in large differences to the single ion properties which de¬ 
termine the crystal electric field and, consequently, the 
anisotropy. 

In the case of some rare earth cations, the magnetic 
ground state is robust under the application of chemical 
pressure. For example, the spin ice state that was first 


observed in Ho 2 Ti 2 Oj 
sure (Ho 2 Sn 2 C >7 24 


rersists at negative chemical pres- 
|) and positive chemical pressure 
(Ho 2 Ge 2 07 [26), [27| ). In other cases, chemical pressure 
can radically alter the magnetic ground state. For exam¬ 
ple, in the ytterbium based pyrochlores, Yb 2 B 2 07 (B = 
Sn, Ti, Ge), the interactions between ytterbium moments 
are very sensitive to changes in chemical pressure [28( . 
Consequently, each of these pyrochlores have a markedly 
different magnetic ground state ranging from ferromag¬ 
netic Coulomb liquid in Yb?Sn?Q 7 [291] to antiferromag¬ 
netic order in Yb 2 Ge 207 [28|. Different magnetic ground 
states are also observed in the terbium based pyrochlores, 
Tb 2 B 2 0 7 (B = Sn, Ti, Ge) @£a,lsl|. 

Even more extreme differences in chemical pressure 
can be achieved through the substitution of lead onto 
the B-site of the pyrochlore lattice. Using high pressure 
synthesis techniques, the pyrochlore stability field was ex¬ 
panded in 1969 by A. W. Sleight to include the plumbates 
1321 . Pb 4+ , with an ionic radius of 77.5 pm, is the largest 
B-site cation to be successfully incorporated into the py¬ 
rochlore structure [32]. However, the very large ionic 
radii of lead imposes limitations on possible A-site rare 
earth cations. Only the largest of the rare earths have 
been successfully prepared: A 2 Pb 2 07 (A = Pr, Nd, Gd). 
In fact, the only rare earth which can be prepared on 
both a titanate and plumbate lattice is gadolinium. In¬ 
stead, the stannate, A 2 Sn 2 07 , and zirconate, A 2 Zr 2 07 , 
pyrochlores, are a more direct point of comparison. Each 
of the above mentioned rare earths can be combined with 
tin and zirconium onto the pyrochlore lattice. 

The praeseodymium pyrochlores, Pr 2 Sn 2 07 and 
Pr 2 Zr 2 07 play host to the so-called dynamic or quan¬ 
tum spin ice states |33[ [34|. Much like the dipolar spin 
ice state, the Pr 3+ moments are locally arranged into 
a two-in/two-out configuration. However, the Pr 3+ mo¬ 
ments are subject to stronger quantum fluctuations and 
remain dynamic down to the lowest measured tempera¬ 
tures, 200 rnK 33]. This state is of significant experimen¬ 


tal 33 h 35| and theoretical interest 36-14lj as it may per¬ 
mit the propagation of monopole excitations. Neutron 
scattering measurements of single crystal Pr 2 Zr 2 07 have 
revealed a large inelastic signal, absent of pinch points, 
suggesting quantum fluctuations of magnetic monopoles 


341. 


There has been comparatively little study of the 
neodymium based pyrochlores. Long range ordering has 
been observed in Nd 2 Sn 2 07 and Nd 2 Zr 2 07 using heat ca¬ 
pacity [H] and magnetic susceptibility [ 20 ]- Nd 2 Sn 2 07 
has a negative Weiss temperature, dew = —0.31 K [20]. 
However, until recently, the ground state of this mate¬ 
rial had not been determined. Recent magnetic neutron 
diffraction measurements have revealed the exact mag¬ 
netic ground state in Nd 2 Sn 2 07 is an antiferromagnetic 
q = 0 state [rj ]. 

Gadolinium pyrochlores are unique in two respects, (i) 
Firstly, gadolinium has a spin only total angular mo¬ 
mentum (S = 7/2, L = 0). As a result, the single 
ion anisotropy of Gd 3+ is small compared to other mag¬ 
netic rare earth cations, (ii) Secondly, gadolinium py¬ 
rochlores have the largest size range for non-magnetic 
B-site cations, ranging from the smallest (Ge 4+ [2l|) to 
the largest (Pb 4+ [32|). However, only the intermedi¬ 
ate sized members, Gd 2 Ti 2 07 and Gd 2 Sn 2 07 , have been 
magnetically characterized until now [l7|, 14414461 ] . Both 
Gd 2 Ti 2 07 and Gd 2 Sn 2 07 are candidates for the Palmer- 
Chalker ground state, a four sublattice Neel state 147 [. 
Gd 2 Sn 2 07 does appear to be a good experimental re¬ 
alization of this state [48]. However, the exact ground 
state of Gd 2 Ti 2 07 has remained experimentally elusive 

[muni. 


With a diverse array of magnetic behavior, the frus¬ 
trated pyrochlores have maintained a high level of inter¬ 
est over the past 20 years. An expanded knowledge of 
their phase diagrams will improve understanding of these 
unique ground states, and may allow the discovery of 
new ground states. Here, we report on the synthesis and 
characterization of the lead-based pyrochlores, A 2 Pb 2 07 
(A = Pr, Nd, Gd). 


EXPERIMENTAL METHODS 

Polycrystalline samples of A 2 Pb 2 07 (A = La, Nd, Pr, 
Gd) were synthesized from stoichiometric amounts of 
A 2 O 3 and Pb0 2 reacted in a gold capsule. The high 
pressure synthesis was carried out in a Walker-type multi- 
anvil press at 2 GPa and 900°C. Room temperature x- 
ray powder diffraction was performed with a D2-phaser 
diffractometer using Cu K-a radiation. Rietveld analy¬ 
ses of the x-ray diffraction patterns were performed using 
Fullprof [5(| . Powder neutron diffraction measurements 
on a 90 mg sample of Pr 2 Pb 2 07 were made on the Gen¬ 
eral Materials Diffractometer (GEM) at the ISIS neutron 
facility. Rietveld analysis of the neutron diffraction pat- 
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terns were performed with GSAS. 

The susceptibility and magnetization measurements 
were made using a Quantum Design MPMS magnetome¬ 
ter equipped with a 3 He insert. The dc magnetic sus¬ 
ceptibility measurements were made in fields of 1000 Oe 
between 0.5 K and 300 K. The ac magnetic suscepti¬ 
bility measurements were carried out with an ac driv¬ 
ing field of 3.5 Oe, over a frequency range of 1 Hz to 
1000 Hz. Specific heat measurements were carried out us¬ 
ing a Quantum Design Physical Property Measurement 
System (PPMS) equipped with a 3 He insert. The mag¬ 
netic component of the specific heat was isolated by sub¬ 
tracting the specific heat of a non-magnetic lattice equiv¬ 
alent, La 2 Pb 207 . 


CRYSTAL STRUCTURES 


TABLE I. Ionic radii for A = Pr, Nd, Gd and B = Ti, Sn, Pb, 
and the resultant pyrochlore lattice parameter. The lattice 
parameters for the stannate pyrochlores are taken from [51] . 


Cation 

Pr (0.99 A) 

Nd (0.98 A) 

Gd (0.94 A) 

Ti (0.61 A) 
Sn (0.69 A) 
Pb (0.78 A) 

10.600 A 
10.872 A 

10.567 A 
10.834 A 

10.18 A 0 - 

10.454 A 
10.729 A 


TABLE II. Selected parameters for the refinement of the neu¬ 
tron diffraction pattern of Pr 2 Pb 2 C >7 


Atom 

Site 

(x,y,z) 

Occupancy 

U iso (AQ 

Pr/Pb 

16d 

(4 A 4) 

V 2 ’ 2 ’ 2 / 

0.924(7)/0.076(7) 

0.0071(9) 

Pb/Pr 

16c 

(0,0,0) 

0.924(7)/0.076(7) 

0.0026(4) 

01 

8 b 

(3 3 3x 

V 8 ’ 8 ’ 8 / 

0.91(1) 

0.006(1) 

02 

48/ 

(0.342,4,4) 

1.0 

0.0085(4) 
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FIG. 1. (Color online) Room temperature powder x-ray 
diffraction patterns and structure refinements of A 2 PI 52 O 7 , 
where A = (a) La, (b) Nd, (c) Pr, and (d) Gd. Patterns for 
the observed (red marks), calculated (black lines), and differ¬ 
ence (blue lines) are shown. The ticks indicate the positions 
of the Bragg reflections for the Fd3m space group. 


The crystal structures of the powder A 2 Pb 2 07 (A = 
La, Nd, Pr, Gd) samples were investigated by room tem¬ 
perature powder x-ray diffraction. The crystal structure 
was determined by Rietveld refinement of the diffraction 
patterns using the space group Fd3m (Figure[l]). Lattice 
parameters for A 2 Pb 2 07 (A = La, Nd, Pr, Gd), listed in 
Table [H were in agreement with prior investigations [32| . 


The lattice parameters for the plumbate pyrochlores fol¬ 
low a linear relation with A site ionic radii. The marked 
peak in each diffraction pattern at 29 = 25.6° corresponds 
to the (222) pyrochlore reflection from the Cu K-/3 wave¬ 
length. 

Lead is the largest known possible B-site cation for 
the pyrochlore lattice. As a result, the A-site to B-site 
ratio, A 3+ /B 4+ , in the lead pyrochlores are amongst the 
smallest. Consequently, much like zirconium-based py¬ 
rochlores, non-negligible site-mix ing between the A and 
B sites is expected to be likely [52]]. To determine the 
exact level of site-mixing in the lead pyrochlores, we 
have performed room temperature neutron diffraction on 
Pr 2 Pb 207 (Figure[2]). Rietveld analysis revealed the pres¬ 
ence of mixing between Pr 3+ and Pb 4+ on the order of 
8 %. This is also evidenced by the reduced size of the 
(111) superlattice peak at 29 = 14° with decreasing rare 
earth size (Figure [1]). The pyrochlore superlattice peaks 
arise due to A and B-site ordering, and vanish in the 
case of random site occupancy. The detailed results of 
the structural refinement of Pi' 2 Pb 207 are presented in 
Table DU 

Gd 2 Pb 207 has an A-site to B-site ratio of 1.36, the 
smallest of any pyrochlore [§]. From inspection of the 
x-ray diffraction pattern, it is apparent that the (111) su¬ 
perlattice peak at 29 = 14° is nearly extinct in Gd 2 Pb 207 
(FigureQJd)). However, neutron diffraction of Gd 2 Pb 207 
is made difficult by the high neutron absorption cross sec¬ 
tion of gadolinium. Instead, we used x-ray refinements to 
investigate the degree of site-mixing in Gd 2 Pb 207 . Re¬ 
finements in which the A-site and B-site occupancies are 
allowed to freely vary give 61% order of the A and B 
sites (39% site mixing), close to the random occupan¬ 
cies (50%) expected for the defect-fluorite structure. For 
comparison, x-ray refinements of the other three lead py¬ 
rochlores we investigated, A 2 Pb 207 (A = La, Nd, Pr) 
each gave 96% order (4% site mixing), which agrees well 
with the neutron refinements. 
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Temperature (K) 


FIG. 3. (Color online) The temperature dependence of the 
magnetic susceptibility (red) and inverse susceptibility (blue) 
of Pr 2 Pb 207 measured with a field of 1000 Oe. Inset: Low 
temperature inverse susceptibility with a Curie-Weiss law fit 
between 0.5 K and 5 K. 


Weiss temperature. The strength of the nearest neighbor 
dipolar interactions, which can be approximated in py- 
rochlores by: 


Arn = 


5/xp n 2 
3 • 47 t r 3 


( 1 ) 


FIG. 2. (Color online) Rietveld refinement of the room tem¬ 
perature TOF diffraction pattern of Pr 2 Pb 2 07 measured on 
the GEM diffractometer. The observed (red marks), calcu¬ 
lated (black lines) and difference (blue lines) patterns are 
shown. The ticks indicate the positions of the Bragg reflec¬ 
tions for the F dim space group. The weighted profile R-factor 
for the fit is R wp = 4.14%. 


MAGNETIC PROPERTIES: Pr 2 Pb 2 0 7 


Figure [3] shows the magnetic susceptibility of 
Pr 2 Pb 207 measured in a field of 1000 Oe. The suscep¬ 
tibility of Pr 2 Pb 2 07 is free of anomalies down to 0.5 K. 
The inverse susceptibility has a large negative curvature 
with increasing temperature, likely due to a large split¬ 
ting of the crystal electric field levels. Consequently, the 
region over which the Curie-Weiss law can be applied 
is limited to low temperatures, between 0.5 K and 5 K 
(Figure [3] inset). The effect of low lying crystal fields 
is to limit the quantitative validity of the Curie-Weiss 
analysis; however, we proceed with the intention of draw¬ 
ing comparisons with other praeseodymium pyrochlores. 
The fit gives a Weiss temperature of dew = —0.74(1) K, 
indicating net antiferromagnetic interactions. The Weiss 
temperatures in Pr 2 Sn 2 07 and Pr 2 Zr 2 07 are 0.32 K [20] 
and —1.41 K 34|, respectively. However, the tempera¬ 
ture regions which exhibit Curie-Weiss behavior differs 
significantly between these three compounds. As a re¬ 
sult, it is difficult to draw comparisons on the basis of the 


when /i is estimated from the low temperature suscepti¬ 
bility, is only 0.11 K in Pr 2 Pb 2 07 . Thus, the magnetism 
in Pr 2 Pb 2 07 is dominated by exchange interactions. 

The higher temperature magnetic susceptibility of 
Pr 2 Pb 2 07 is dominated by single ion physics, due to a 
large crystal electric field splitting. The 9-fold degener¬ 
ate state of Pr 3+ is split into three non-Kramers doublets 
and three singlets [20|. From the Curie-Weiss fit between 
0.5 and 5 K, we calculate an effective magnetic moment 
of 2.53(1) /ie, a value not significantly different from that 
found in Pr 2 Sn 2 07 and Pr 2 Zr 2 07 [20|,|34j]. However, this 
moment falls short of the 3.58 /ie free-ion moment for a 
pure |J = ±4) state, suggesting a mixing of |J = ±1) 
and |J = =f 2) terms j2c|. Thus, Pr 2 Pb 2 07 has a doublet 
ground state with Ising moments pointing along the local 
< 111 > trigonal axes. 

The magnetization of Pr 2 Pb 2 07 , shown in Figure [4j 
does not saturate in a field of 7 T at 0.5 K. This behavior 
deviates strongly from the Brillouin function for J = 4, 

. . 2 J + 1 ( 2 J + 1 \ 1 / x \ . . 

Bj(l) = ~ 2 T\~ZTV - Tj coUl (27) (2) 

where rr = (3) 

k B T 

due to significant interactions between the Pr 3+ moments 
at these temperatures and deviations from isotropic crys¬ 
tal field levels. The magnetization of Pr 2 Pb 2 C >7 appears 
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FIG. 4. (Color online) The magnetization of Pr 2 Pb 2 07 as a 
function of applied field between 0 T and 7 T at 0.5 K and 
1.8 K. The dashed line marks half the effective moment value, 
1.25 /r B /Pr. 



Temperature (K) 

FIG. 5. (Color online) The specific heat of Pr 2 Pb 2 C >7 and 
non-magnetic La 2 Pb 207 . Inset: low temperature lattice sub¬ 
tracted specific heat. 


to be approaching a saturation value of ~1.25 hb = 

A saturation value of half the effective moment in pow¬ 
der samples is a hallmark of < 111 > Ising anisot ropy , as 
observed in other praseodymium pyrochlores |53l l54| . In 
general, the magnetization of Pr 2 Pb 2 C >7 closely resem¬ 
bles that of Pr 2 Zr 2 C> 7 , which also fails to saturate by 7 T 
541. 

The heat capacity of Pr 2 Pb 207 contains a broad 
anomaly centered at 1.2 K (Figure [5]). This feature does 
not resemble the sharp anomalies which typically signify 
long-range ordering, nor does it fit to a Schottky model. 
Instead, the peak in the heat capacity of Pr 2 Pb 2 07 bears 
a resemblance to those observed in spin ice, attributed 
to a condensation of thermally activated monopole pairs 
Eel . Much like Pr 2 Zr 207 , the peak is broader than those 
observed in the dipolar spin ices [i]. This broadening is 
attributed to monopole quantum dynamics in Pr 2 Zr 207 , 
and is likely also at play in Pr 2 Pb 2 C> 7 . 

The lattice contribution to the specific heat was 
subtracted by using the specific heat of non-magnetic 
La 2 Pb 2 C> 7 . Calculation of the entropy for the lattice sub¬ 
tracted specific heat, AS = f dT , gives a value which 
exceeds the maximum value for an Ising system, Rln( 2). 
However, there are other known contributions to the spe¬ 
cific heat of Pr 2 Pb 2 C> 7 . The hyperfine splitting of the nu¬ 
clear levels of Pr 3+ produces a nuclear Schottky anomaly 
with significant intensity below 1 K. Furthermore, the 
crystal field scheme of Pr 2 Pb 2 07 has not been measured, 
but there is likely a low-lying crystal field contributing to 
the specific heat, as evidenced from changes in the curva¬ 
ture close to 10 K. However, to subtract this contribution 
would require a more precise knowledge of the low-lying 
crystal field levels. Since these contributions to the spe¬ 
cific heat of Pr 2 Pb 2 C >7 cannot be accurately determined, 
a realistic calculation of the magnetic entropy is not pos¬ 


sible. It is interesting to note that the lowest crystal fields 
in Pi^S^CH and Pr^Zi^CH occur at energies an order of 
magnitude higher [34|, |55|. Thus, Pr 2 Pb 207 represents 
an interesting opportunity to investigate the impact of a 
low-lying crystal field on a dynamic spin ice state. 


MAGNETIC PROPERTIES: Nd 2 Pb 2 0 7 

The magnetic susceptibility of Nd 2 Pb 2 C> 7 , shown in 
Figure 0 contains no anomalies down to 0.5 K. Between 
0.5 K and 15 K, the inverse susceptibility is well-fit by the 
Curie-Weiss law. Before proceeding, we should note that 
much like Pr 2 Pb 2 C> 7 , the crystal electric field is likely im¬ 
pacting the analysis. The Weiss temperature calculated 
from this fit is 0cw = —0.069(4) K, indicating weak net 
antiferromagnetic interactions. A Curie-Weiss fit over 
the same temperature range in Nd 2 Sn 2 C >7 gives a Weiss 
temperature of —0.31 K [20j. Thus, the change from the 
stannate to the plumbate results in a very close competi¬ 
tion of antiferromagnetic and ferromagnetic interactions. 
It is worth noting that the calculated dipolar interactions 
in Nd 2 Pb 2 C >7 and INchS^CH are 0.11 K and 0.12 K, re¬ 
spectively. Thus, the changing strength of the dipolar 
interaction alone cannot account for the change in the 
Weiss temperature. The decreased Weiss temperature in 
Nd 2 Pb 2 07 is thus indicative of decreased exchange inter¬ 
actions. 

The calculated effective magnetic moment for 
Nd 2 Pb 2 07 , from the Curie-Weiss fit, is 2.55(7) ^b- 
This value falls short of the free-ion value of 3.62 
for | J = ±9/2), indicating the mixing of | J = ±3/2) 
terms. Above 20 K, the susceptibility deviates from 
Curie-Weiss behavior due to low lying crystal field levels. 
The 10-fold degenerate ground state of Nd 3+ (J = 9/2) 
is split into five symmetry protected Kramers doublets. 
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Temperature (K) 

FIG. 6. (Color online) The magnetic susceptibility of 
Nd 2 Pb 2 0 7 measured with an applied field of 1000 Oe. The 
left axis gives the direct susceptibility (red) and the right axis 
gives the inverse susceptibility (blue). The Curie-Weiss fit 
between 0.5 K and 15 K is given by the blue line. 



FIG. 8. (Color online) Magnetic specific heat of Nd 2 Pb 2 0 7 
measured in zero field, 3 T, and 9 T between 0.4 K and 40 K. 
Inset: Scaled in view of the low temperature zero field specific 
heat of Nd 2 Pb 2 0 7 . 



FIG. 7. (Color online) The magnetization of Nd 2 Pb 2 0 7 as a 
function of applied field between 0 T and 7 T at 0.5 K and 
2 K. The solid blue and red curves give the Brillouin function 
for 0.5 K and 2 K respectively. 



Temperature (K) 


FIG. 9. (Color online) Calculated magnetic entropy of 
Nd 2 Pb 2 0 7 measured in zero field, 3 T, and 9 T between 
0.4 K and 40 K. The dashed line indicates the expected value, 
I? In (2), for an Ising system. 


Figure [3 shows the magnetization of Nd2Pb20 7 at 0.5 K 
and 2 K. At both temperatures, the magnetization is 
similar to that of a Brillouin function for J = 9/2. The 
magnetization saturates at 1.26 /XB/Nd 3+ , half the value 
of the effective moment. Again, this is a consequence of 
local < 111 > Ising anisotropy and powder averaging. 

The magnetic specific heat of Nd2Pb20 7 is shown in 
Figure [8] In zero field, the magnetic component of the 
specific heat increases below 3 K. Below 0.5 K, the spe¬ 
cific heat begins to strongly increase, with its maximum 
occurring at 0.41 K. However, only a single data point 
could be collected below the peak’s maximum. This peak, 
which is likely due to long-range magnetic ordering, shifts 
upwards in temperature and broadens with an externally 
applied field. A field of 3 T shifts the maximum to 2.5 K, 


while a field of 9 T further increases it to 8 K. The mag¬ 
netic origin of this peak is most easily verified by inte¬ 
grating out the entropy, which should give R ■ In 2 for 
an Ising system. However, this cannot be done with our 
zero-field measurement, in which the decreasing side of 
the anomaly is outside our measurement range (Figure [SJ. 
Integration of C p /T between 0.4 K and 40 K for the 3 T 
and 9 T data sets yields 0.95-1? In 2 and 0.96-1?In 2, re¬ 
spectively (Figure [9]). 

Long-range order was also found in Nd2Sii20 7 , which 
has both a sharp peak in its heat capacity [42| and a cusp 
in its susceptibility [2Cj at 0.9 K. It is interesting to note 
that the ordering temperature in Nd 2 Pb 2 0 7 is reduced 
from Nd 2 Sn 2 0 7 by a factor of two. Recent neutron scat¬ 
tering measurements of Nd 2 Sn 20 7 have confirmed that 
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the ordered state is the k=0 antiferromagnetic all-in/all- 
out state [43]. Neutron scattering will also be required 
to determine the consequence of net ferromagnetic inter¬ 
actions on the magnetic ground state of Nd 2 Pb 2 0 7 . 


MAGNETIC PROPERTIES: Gd 2 Pb 2 0 7 


Figure HUl' ai shows the low temperature magnetic sus¬ 
ceptibility of Gd 2 Pb 2 0 7 measured in 1000 Oe. A mag¬ 
netic ordering transition is signaled by a cusp in the 
susceptibility at 0.8 K. There is a slight splitting of the 
field cooled and zero field cooled susceptibility below the 
transition at 0.8 K indicating a small amount of irre¬ 
versibility. The susceptibilities of the related pyrochlores, 
Gd 2 Ti 2 0 7 and Gd 2 Sn 2 0 7 , both contain magnetic order¬ 
ing anomalies close to 1.0 K 0 - They too both exhibit 
a degree of irreversibility below these transitions. The 
maximum FC/ZFC splitting measured in Gd 2 Pb 2 0 7 is 
0.03 emu/mol-Gd 3+ , v ery similar to the magnitude of the 
splitting in Gd 2 Ti 2 0 7 [4J]. Conversely, in Gd 2 Sn 2 0 7 , the 
FC/ZFC splitting is an order of magnitude larger 44]. 

The high temperature magnetic susceptibility of 
Gd 2 Pb 2 0 7 up to 300 K is shown in Figure [ToT bl. Be¬ 
tween 10 K and 300 K the inverse susceptibility is strik¬ 
ingly linear, fitting well to Curie-Weiss behavior. The 
Weiss temperature is —7.38(6) K, indicating dominant an¬ 
tiferromagnetic interactions. The Weiss temperatures for 
Gd 2 Ti 2 0 7 and Cd 2 Sn 2 0 7 are —9.9(1) K and —8.6(1) K, 
respectively 44]. As anticipated, the substitution of the 
larger Pb 4+ cation appears to have weakened the antifer¬ 
romagnetic exchange in Gd 2 Pb 2 0 7 . 

The Curie-Weiss fit of Gd 2 Pb 2 0 7 can also be used to 
extract an effective moment, yielding 7.67(1) hb (Figure 
M b))- This calculated moment is close to the expected 
7.94 /i b free-ion moment. In Gd 3+ , the isotropic spin, 
S = 7/2, is the only contribution to the total angu¬ 
lar momentum, J. Consequently, the 8-fold degeneracy 
of the ground state is only marginally lifted by the crys¬ 
tal electric field, on the order of a fraction of a Kelvin 


17]. The high degree of linearity in the inverse suscepti¬ 


bility between 10 K to 300 K suggests there are no well- 
separated crystal field levels in Gd 2 Pb 2 0 7 , supporting 
the picture described above. 

The magnetization of Gd 2 Pb 2 0 7 is shown in FigureiJTl 
at temperatures ranging from 0.5 K to 50 K. At 10 K 
and higher, the magnetization fits well to the Brillouin 
function for J = 7/2, which describes a system of para¬ 
magnetic, non-interacting moments. At 0.5 K, where the 
Brillouin function fails completely, the magnetization ap¬ 
proaches a saturation value on the order of the free-ion 
Gd 3+ moment by 7 T. 

The magnetic specific heat of Gd 2 Pb 2 0 7 is shown in 
Figure [T2] The magnetic specific heat increases below 
10 K and peaks at 0.81 K. Extrapolating C mag /T to 
0 J/mol-K 2 , gives an integrated entropy of 17.1 J/mol- 




FIG. 10. (Color online) (a) The low temperature magnetic 
susceptibility of Gd 2 Pb 2 0 7 measured with an applied field of 
1000 Oe. Below the magnetic ordering transition at 0.8 K 
there is a splitting between the zero field cooled and field 
cooled susceptibility, (b) The high temperature magnetic sus¬ 
ceptibility of Gd 2 Pb20 7 . The left axis gives the direct suscep¬ 
tibility (red) and the right axis gives the inverse susceptibility 
(blue). The Curie-Weiss law fit is given by the blue line. 



FIG. 11. (Color online) The magnetization of Gd 2 Pb 2 0 7 as a 
function of applied field between 0 T and 7 T at temperatures 
between 0.5 K and 50 K. The solid curves show the Brillouin 
function for J = 7/2 at each corresponding temperature. 
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Temperature (K) 

FIG. 12. (Color online) Magnetic specific heat of Gd 2 Pb 207 
measured in zero field between 0.4 K and 20 K. Inset: In¬ 
tegrated entropy of Gd 2 Pb 207 , giving the expected value, 
I?. In 8, for Heisenberg-like spins. 



Temperature (K) 

FIG. 13. (Color online) Power law fits to the low temperature 
specific heat of Gd 2 Pb 207 . The least squares refinement gives 
a Cmag oc T 3,/2 dependence (red curve). 


K (Figure [T3] inset). This value corresponds to 
0.99-i?ln8 = 0.99-i?ln (25 + 1). This is the expected 
value for Gd 3+ (S = 7/2), which, with very small single 
ion anisotropy, has Heisenberg-like spins. Anomalies re¬ 
lated to magnetic ordering are also observed in Gd 2 Ti 2 C >7 
and Gd 2 Sn 207 at approximately 1 K. The maximum am¬ 
plitude of the anomaly in Gd 2 Pb 207 is 8.7 J/mol-K. This 
is similar in magnitude to that of Gd 2 Ti 2 07 , whereas 
the anomaly in Gd 2 Sn 2 07 has a maximum amplitude 
of 120 J/mol-K [Sj. However, in Gd 2 Ti 207 , there is a 
second anomaly at 0.75 K and the order ing in this com¬ 
pound occurs in two distinct stages 44j, |4(j. Thus, the 
heat capacity anomaly of Gd 2 Pb 207 does not strongly 
resemble either of its group members. 

The temperature dependence of the magnetic specific 
heat of Gd 2 Pb 207 , below the ordering temperature, re- 


TABLE III. Comparison of lattice parameters and selected 
magnetic parameters in the A 2 Pb 207 (A = La, Nd, Pr and 
Gd pyrochlores). 



a (A) 

#cw (K) 

Meff (Mb) 

Anisotropy 

La 2 Pb 2 07 

10.9682(4) 

- 

- 

- 

Pr2Pb207 

10.8721(9) 

—0.74(1) 

2.53(1) 

Ising 

Nd2Pb207 

10.8336(4) 

0.15(1) 

2.47(1) 

Ising 

Gd2Pb207 

10.7292(8) 

-7.38(6) 

7.67(1) 

Heisenberg 


veals further differences from the other gadolinium py¬ 
rochlores (Figure fl3l) . Conventional ungapped antiferro¬ 
magnetic spin wave excitations give rise to a C m&g oc T 3 . 
However, the specific heat below the ordering transition 
in both Gd 2 Ti 2 C >7 and Gd 2 Sn 2 C >7 is better parameterized 


by Gmag oc T 2 [4J,|56|. In Gd 2 Sn 207 , the specific heat is 
found to decrease exponentially below 350 mK 57]. This 
behavior in Gd 2 Sn 2 07 is attributed to conventional anti¬ 
ferromagnetic gapped magnons and can be theoretically 
reconciled with linear spin-wave theory (57| . However, 
the Cmag oc T 2 dependence in Gd 2 Ti 207 , which is un¬ 
gapped down to 100 mK, remains a mystery [561. 

A least squares Et to the heat capacity of Gd 2 Pb 2 C >7 
between 0.35 K and 0.75 K to G mag oc T v gives 
?7 = 1.49(3). Thus, the specific heat of Gd 2 Pb 2 C >7 be¬ 
low 0.8 K is well described by a G mag oc T 3 / 2 power 
law, as shown in Figure 1131 In an ordered ferromag- 
net, spin wave excitations make a C ma g oc T 3 / 2 con¬ 
tribution to the specific heat. However, a ferromagnetic 
ground state cannot be reconciled with the overwhelm¬ 
ingly negative (antiferromagnetic) Curie-Weiss temper¬ 
ature in Gd 2 Pb 207 . Determination of the magnetic 
ground state of Gd 2 Pb 2 C >7 will require neutron diffrac¬ 
tion, which is made challenging by the high neutron 
absorption cross section of gadolinium. While the ex¬ 
act ground state remains unknown at present, it is clear 
that it is markedly different from those of Gd 2 Ti 207 and 
Gd 2 Sn 207 . In the absence of neutron diffraction or muon 
spin relaxation microscopic magnetic measurements we 
cannot rule out spin glass freezing. However, the sharp¬ 
ness of the transition seen in our bulk dc suscepbitility 
and specific heat measurements, along with a lack of fre¬ 
quency dependence in ac susceptibility (not shown) make 
spin glass ordering highly unlikely in Gd 2 Pb 2 C> 7 . This is 
particularily noteworthy given the significant degree of 
site mixing. 


CONCLUSIONS 

We have reported on the crystallographic and mag¬ 
netic properties of the lead based pyrochlores A 2 Pb 2 C >7 
(A = La, Pr, Nd, Gd) (Table Mil) . We have explored 
a new region of phase space in a family of materials 
known to have a rich collection of ground states. Lead 
is the largest B-site cation which can occupy the py- 
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rochlore lattice. Consequently, lead-substitution on the 
pyrochlore lattice results in an effective weakening of the 
exchange interaction. As a result, magnetic ordering is 
suppressed to lower temperatures in each case. Pr 2 Pb 2 07 
has strong Ising < 111 > anisotropy and an absence of 
long-range magnetic order down to 500 mK. Pr 2 Pb 2 07 
has a spin ice-like anomaly in its heat capacity centered 
at 1.2 K. Thus, Pr 2 Pb 2 07 is a candidate for a quantum 
spin ice state, despite weakened exchange and mild disor¬ 
der. Nd 2 Pb 2 07 also has strong <111> Ising anisotropy. 
However, Nd 2 Pb 2 07 transitions to long-range magnetic 
order at 0.41 K. The small, negative Weiss temperature 
of Nd 2 Pb 2 07 indicates it may lie on the border between 
ferromagnetism and antiferromagnetism. Gd 2 Pb 2 07 is 
a good approximation to Heisenberg spins on the py¬ 
rochlore lattice. However, due to the very similar ionic 
radii of gadolinium and lead, Gd 2 Pb 2 07 has a significant 
level of site mixing. Despite this structural disorder, mag¬ 
netic ordering in Gd 2 Pb 2 07 at 0.81 K is signaled by a 
cusp in the susceptibility and a lambda-like specific heat 
anomaly. Unlike Gd 2 Ti 2 C >7 the ordering in Gd 2 Pb 2 07 
occurs in a single stage. Furthermore, a T 3 1 4 5 / 2 depen¬ 
dence is observed in the low temperature specific heat, 
a value which is typically associated with the spin wave 
contribution in a conventional ferromagnet. Thus, it is 
unlikely that Gd 2 Pb 2 07 is adopting the antiferromag¬ 
netic Palmer-Chalker ground state. To that end, neu¬ 
tron diffraction of Gd 2 Pb 2 C >7 is recommended. In gen¬ 
eral, these lead-based pyrochlores provide insight into the 
effects of weaker exchange on highly frustrated ground 
states. 
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